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Several Wilson loops on several lattice sizes are computed in Perturbation Theory via a stochastic method. 
Applications include: Renormalons, the Mass Term in HQET and (possibly) the /3-function. 



Outlook 

Wilson Loops (WL) were the historic play- 
ground (and success . . . ) of Numerical Stochastic 
Perturbation Theory (NSPT) for Lattice Gauge 
Theory jy. Having by now an increased com- 
puting power available, we are computing high 
perturbative orders of various WL on various lat- 
tice sizes in Lattice SU(3). Physical motiva- 
tions range over a variety of issues (not every one 
within reach, at the moment): Renormalons and 
Lattice Perturbation Theory (LPT), the Mass 
Term in Heavy Quark Effective Theory (HQET) 
and the Lattice /3-function. 

1. Renormalons and LPT 

In g| WL of sizes lxl and 2x2 were computed 
in LPT via NSPT up to /3" 8 order. The expected 
Rcnormalon contribution was found according to 
the formula (Q = a -1 ) 
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which is fixed by dimensional and Renormaliza- 
tion Group considerations. (0) can be cast in a 
form from which a power expansion can be eas- 
ily extracted (two loop asymptotic scaling is as- 
sumed) 
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Actually (Q) refers to some continuum scheme. 
On a finite lattice (which is what NSPT needs) 
one has to deal with 
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The factor s in the argument of a s is in charge of 
the lattice-continuum matching and is expected 

to be of the order s ~ (^ \ lat t ) with respect to 
some continuum scheme. An explicit IR cut-off is 
present, dependent on the lattice size (Qo(N) = 
2ir{Na)~ 1 , where N is the number of points in 
any direction). (|^) results in a new power expan- 
sion 
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whose coefficients C^ on (C, s, N) are given in terms 
of incomplete T-functions. In 0] (via a slightly 
different, equivalent formalism) it was shown that 
this renormalon contribution can account for the 
growth of the first 8 coefficients in the pertubative 
expansion of the plaquette 
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1. e. convenient values C* and s* were fitted so 
that c] xl were recognized to be asymptotically 
the same as C™ n (C*, s*, 8) (in (2| computations 
were performed on a rather small N — 8 lat- 
tice); Cf n (C**,s*,8) (with a different choice for 
C) were also shown to fit the expansion of W 2x2 . 

Control on this (renormalon) perturbative con- 
tribution is crucial in the analysis presented in || : 
once its resummation is subtracted from Monte 
Carlo data for the plaquette, one is left with a 
0(A 2 /Q 2 ) contribution rather than the expected 
0(A 4 /Q 4 ). In order to trust this remarkable re- 
sult, one would like to further test the asymptotic 
formula ([|) by going to even higher orders in the 
perturbative expansion on various lattice sizes. 

By now we know the expansion of the basic 
plaquette up to order fi~ w on both N = 8 and 
N = 24 lattices. We do not present here the 
definitive results (which will be published soon 
elsewhere), but we do show how well the new re- 
sults are described by C; on (C*, a*, N) {i.e. by the 
values for C and s that were fitted in ||||). I n the 
figure the expected Cf n (C*, a*, 24) are plotted 
together with the computed c] xl (£ — 1 . . . 10) 
for the N = 24 lattice. The renormalon contri- 
bution is indeed there just like described in 0. 
Work is in progress to perform the whole analysis 
of |] on the N = 24 lattice. 

2. The Mass Term in HQET 

Consider WL of various sizes (in particular 
square L x L loops). From their renormalization 
properties Q one has to expect 

W(L) =cxp- 4Mi w(L) (4) 

The first factor is the exponential of the perime- 
ter times the Mass Term one has to deal with 
in HQET (an additive, linearly divergent mass 
renormalization). w(L) contains logarithmic di- 
vergences, in particular only those connected to 
the coupling renormalization is there is no "cor- 
ner" (which is not the case on the lattice). As 
stressed in Hashimoto's review at this conference 
H (after |J), the mass term is a fundamental 
building block in a renormalon-safe determina- 
tion of the b-quark mass from the lattice. 
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Figure 1. c] xl obtained on a N = 24 lattice (cir- 
cles) and C r e cn (C*,s*,24) (crosses) vs perturba- 
tive order I. 



The mass term could of course be determined 
from the heavy quark propagator. The determi- 
nation we are going to report on (again after f|) 
goes through the computation of various WL and 
is well suited for NSPT as it is gauge invariant Q]. 
By computing perturbative expansions of WL of 
various sizes 

W(L) = 1 - Wi(L) /T 1 - W 2 (L) [3~ 2 + ... 
and therefore 

-\ogW{L) = Wi(L)p- 1 + 

(V 2 (L) + iwMi) 2 ) r 2 + ... 

one can extract from each order the leading (lin- 
ear) behaviour in L, that is mass term 

4 M L = 4 (Mi /T 1 + M 2 /T 2 + . . .) L . (5) 

At the moment we have got results for square 
loops on a (N = 24) 4 lattice. Results on bigger 
lattices (N — 28, 32) and rectangular loops are 
expected to come soon. They will be crucial to 
control finite size effects and subleading (logarith- 
mic) contributions. As Mi and M 2 are already 
known, they have to be recovered. Consider for 
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example Mi , which has to be recovered by fitting 
W X (L) to 

Wi{L) — 4 Mi L + k\ logL + fc 

The analytical result is Mi = 1.01 and we get 
Mi — 1.00 ± 0.02. From |j] one learns the second 
coefficient M% = 2.54, while from our fits we get 
M 2 = 2.43 ± 0.15. The errors we quote depend 
on both finite size effects and fitting subleading 
contributions. Note that since Ws(L) contains 
a (logL) 2 contribution, the latter point is even 
more crucial in the determination of Ms- At the 
moment we can pin down a preliminary number 
M3 = 8±1 which will turn in a definite number as 
soon as we get definite results not only for square 
loops, but also for ratios of rectangular loops. 

3. The Lattice /3-function 

We now turn to describe what could possibly 
be another application of our perturbative com- 
putations. Ratios of WL combined in such a way 
that the corner and mass contributions cancel out 
were introduced several years ago by Creutz || 
in order to study (among other things) the non- 
perturbative Lattice /3-function. 

Consider rectangular WL of sizes (L,L), 
(2L,L), (2L,2L) on a (4L) 4 lattice and form the 
ratio 

W(L, L) W(2L, 2L) 



K(L) = 



W{L,2L) 2 
1 + ci(L) a Q + c 2 (L) a 2 + 



where ao is the bare lattice coupling. The only 
scale is L and the only renormalization needed is 
that of the coupling, so that (for example) 

a n (L) = -±—\og{K{L)) (6) 
ci{L) 

is supposed to be an acceptable coupling running 
with L. Since the pertubative matching between 
coupling in different schemes 

ot\{sn) = a 2 (^) + di(s)a 2 {n) 2 + ... 

contains all the informations about the pertur- 
bative /3-functions in both schemes, one could 
in principle study the perturbative Lattice (3- 
function by computing the perturbative expan- 
sion of acn(L) in a at different values of N, 



where L = Na (i.e. by computing different WL 
on different lattice sizes). 

Within this application the big issue is accu- 
racy. Due to large cancellations of the mass con- 
tribution, a fraction of per mille accuracy on the 
Ci(L) (which is within reach) straight away de- 
generates when one computes for example di(s) 
(still we were able to recover bo, the first universal 
coefficient of the /3-function) . In view of this, it 
is quite unlikely that one can attain the terrific 
accuracy which should be needed to get the first 
unknown information. One should most probably 
look for a smarter definition of the coupling. 

Conclusions 

Preliminary results were reported, coming from 
the computations of perturbative expansions of 
various WL. The basic plaquette is now known 
up to (3~ w order and the results on Renormalons 
^ are indeed to be trusted. Work is in progress 
to gain further confidence in the results of |3j as 
well. A preliminary results on third order in the 
computation of the mass term has been reported 
and a definite result will be published soon. In 
principle also the perturbative lattice /3-function 
could be studied via WL in NSPT, even if the 
accuracy needed to get the first unknown result 
is quite unlikely to be attained. 
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